SLAC-PUB-12867, ZU-TH 25/07 



S@M, a Mathematica Implementation of the 
Spinor-Helicity Formalism 

D. Maitre^'^ , P. Mastrolia^ 

^Institut fiir Theoretische Physik 
University of Zurich 
Winterthurerstrasse 190, CH-8057 Zurich 

^Stanford Linear Accelerator Center 
2575 Sand Hill Road 
Menlo Park, CA 94025 



Abstract 

In this paper we present the package S@M (Spinors@Mathematica) which imple- 
ments the spinor-hehcity formahsm in Mathematica. The package ahows the use of 
complex-spinor algebra along with the multi-purpose features of Mathematica. The 
package defines the spinor objects with their basic properties along with functions 
to manipulate them. It also offers the possibility of evaluating the spinorial objects 
numerically at every computational step. The package is therefore well suited to 
be used in the context of on-shell technology, in particular for the evaluation of 
scattering amplitudes at tree- and loop-level. 
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LONG WRITE-UP 



1 Introduction 

Theoretical understanding of background processes is essential to single out 
interesting signals in the rich landscape of events which will take place at 
the forthcoming CERN experiment, the Large Hadron Collider (LHC). Many 
methods are available for computing Standard Model backgrounds at the lead- 
ing order (LO) in perturbation theory [1-5] , based either on automatic sum- 
mation of tree- level Fcynman diagram, or off-shell recursive algorithms for 
currents [6-8]. But quantitative estimates for most processes require a calcu- 
lation with next-to-leading order (NLO) accuracy - sec for instance [9]. NLO 
calculations require knowledge of both virtual and real-radiation corrections 
to the LO process. While the real-radiation corrections can be computed us- 
ing tree-level techniques, the bottleneck for the availability of results with 
NLO level accuracy [10-19] is the non-trivial evaluation of one-loop virtual 
corrections. New approaches tackling the evaluation of one-loop multi-parton 
amplitudes have recently been under intense development [16, 17,20-51]. 

The spinor-helicity formalism [52] for scattering amplitudes has proven an in- 
valuable tool in perturbative computation since its development in the 1980's, 
being responsible for the discovery of compact representations of tree and 
loop amplitudes. Instead of Lorentz inner products of momenta, it relies on 
the more fundamental spinor products. These neatly capture the analytic 
properties of on-shell scattering amplitudes, like the factorization behavior 
on multi-particle-channels. The recent boost in the progress of evaluating on- 
shell scattering amplitudes is due to turning qualitative information on their 
analytic properties into quantitative tools for computing them. 

On-shell methods [53] restrict the propagating states to the physical ones and 
the spinor-helicity formalism is therefore well suited to avoid the (intermedi- 
ate) treatment of unphysical degrees of freedom whose effects disappear from 
final results. Moreover, on-shell methods are tailored for the parallel treat- 
ment of sets of diagrams which share a common kinematic structure, such as 
multi-particle poles at tree-level and branch-cuts at loop-level [54,55]. They 
arc therefore suitable for extracting analytic information from simpler am- 
plitudes in a recursive/iterative fashion, since the singularities of scattering 
amplitude are determined by lower-point amplitudes in the case of poles and 
by lower- loop ones in the case of cuts [56-58] . 

On-shell methods were originally used in [59] and in the more recent sys- 
tematized implementations for the completion of all six-gluon helicity am- 
plitudes [30,35,36,39,41,42,46-48] and the calculation of the six-photon 
amplitudes [31, 60] in agreement with the numerical results of [29, 33] and 
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[61,62] respectively. 



The singularity information can be extracted by defining amplitudes for suit- 
able complex, yet on-shell, values of the external momenta - an idea that 
stemmed from Witten's development of twistor string theory [63-66] . Generat- 
ing complex momenta by modifying spinor variables, considered as fundamen- 
tal objects, leads to new ways to exploit the kinematic properties of helicity 
amplitudes. The new-born complex momenta have the property of preserving 
overall momentum conservation and on-shell nature. Complex kinematics al- 
low the exploration of singularities of on-shell amplitudes and the use of factor- 
ization information to reconstruct tree amplitudes recursively from their poles. 
The application of factorization in the on-shell method is realized through glu- 
ing lower-point tree amplitudes to form higher points ones linked by on-shell 
yet complex propagating particles. The construction of tree amplitudes via 
on-shell recursion essentially amounts to a reversal of the collinear limit. This 
is made possible by complexifying momenta in their spinorial representation, 
and results in a quadratic recursion, the BCFW recurrence relation [67], which 
works for massive particles [68-71] as well. 

Complex kinematics are useful for the fulfillment of generalized unitarity con- 
ditions as well. At one loop, generalized unitarity corresponds to requiring 
more than two internal particles to be on-shell, such constraints cannot be 
realized in general with real Minkowski momenta. 

The application of unitarity as an on-shell method of calculation is based 
on two principles: i) sewing tree amplitudes together to form one- loop am- 
plitudes; ii) decomposing loop-amplitudes in terms of a basis of scalar loop- 
integrals [72,73]. Matching the generalized cuts of the amplitude with the cuts 
of basic integrals provides an efficient way to extract the rational coefficients 
from the decomposition. The unitarity method [35, 36] provides a technique 
for producing functions with the correct branch cuts in all channels [74], as 
determined by products of tree amplitudes. 

The use of four-dimensional states and momenta in the cuts enable the con- 
struction of the (poly) logarithmic terms in the amplitudes, but generically 
drops rational terms, which have to be recovered independently. 
More recent improvements to the unitarity method [40] use complex momenta 
within generalized unitarity, allowing for a simple and purely algebraic deter- 
mination of box integral coefficients from quadruple-cuts. Using double- and 
triple-unit arity cuts have led to very efficient techniques for extracting trian- 
gle and bubble integral coefficients analytically [42,48,50,75]. In particular 
in [42,48,50] the phase-space integration has been reduced to the extraction 
of residues in spinor variables and, at the occurrence, to trivial Feynman- 
parametric integration. This approach has been used to compute analytically 
some contributions to the six-gluon amplitude [42,48], and the calculation 
of the complete six- photon amplitudes [60], whose cut-constructibility was 
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shown in [31] . Other approaches have combined the knowledge of the generic 
structure of loop-integrand [25, 76] with the simphfication induced by cut- 
constraints, ending up with a unitarity-motivated loop-integral decomposition 
[32,33,62,75]. 

The four-dimensional version of the unitarity method leaves the pure rational- 
function terms in the amplitudes undetermined. New approaches to computing 
rational terms use an optimized organization of Feynman diagrams, by focus- 
ing the standard tensor reduction to tensor integrals which could generate the 
rational terms [30,31]. 

A recent investigation [77] on the source of rational terms has been exploring 
the idea of their generation via a set of Lorentz- violating counterterms. 
Alternatively, these rational functions can be characterized by their kinematic 
poles. An efficient means for constructing these terms from their poles and 
residues is based on BCFW-hke recursion relations [43-47]. 
The rational parts of amplitudes can be also detected with D-dimensional 
unitarity cuts [37,49,58,78,79], and in [49,51] the benefits of four- dimensional 
spinorial integration [42,48,50] have been extended to work within the dimen- 
sional regularization scheme and with massive particles. 

In this paper we present the package S@M (Spinors@iMathematica) which 
implements the spinor-helicity formalism in Mathematica. The package allows 
the use of complex-spinor algebra along with the multi-purpose features of 
Mathematica. The package provides 

• the definitions of the spinor objects with their basic properties, 

• functions to manipulate them 

• numerical evaluation. 

These capabilities make the package S@M suitable for, for example, 

- the generation of complex spinors associated with solutions of multi-particle 
factorization and of generalized-cut conditions; 

- the implementation of BCFW-like recurrence relations for constructing high- 
multiphcity tree amplitudes [67,68] and rational coefficients [43-45,80]; 

- the decomposition of massive momenta onto massless ones, useful for the 
implementation of the MHV- rules [81,82] and for the BCFW-shift of mas- 
sive legs [70]; 

- the algebraic manipulation of products of tree amplitudes with complex 
spinors sewn in unitarity-cuts. 

The tool presented here is therefore oriented towards the evaluation of hehcity 
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amplitudes at LO and beyond, which are relevant for the phenomenology of 
the Standard Model, for the study of the so called constructible theories [83], 
and for the investigation of the ultraviolet-behavior of gravity (see [84, 85] 
and references therein) . 

The paper is organized as follows. Section 2 defines the notation used by the 
package. Section 3 describes the implementation. In last section, we show three 
more involved examples of the use of the package SOM. 



2 Notation 



This section describes the conventions and notation used in the package. 



2. 1 Two-dimensional Spinors 



The two dimensional spinors A and A for a massless fermion with four mo- 
mentum p are defined through the Dirac equation 



^Xip) = , ~X{p)^ = , 



with 



(1) 



1 

1 



a' 



^ p_ —p^ 

^-pi p+ 

-i 

1 



a 



1 
-1 



(2) 
(3) 



The solutions of opposite chirality are 



K{p) = c 



(p. 

[pi 



j and Xa(p) = c (p+ p^^ . 



(4) 



We define A and A as the spinor with up-indices, obtained by contracting the 
spinors A and A with the e-tensor. 
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pz. 



(5) 



where e"'' — e""^ — ia2- 

The underhned notation for spinors carrying up-indices has been introduced 
here to avoid the introduction of indices in the package S@M. By making the 
normahzation choice 



1 



c — c — 



(6) 



we have 



= A,A„ , = A"A«, (7) 

so that the expressions for the two-dimensional spinors read, 



Aa(p) = 



and Aa(p)= 



P- 



(8) 



The the above formulas hold as well for spinors associated to any massless 
complex four- vectors. For massless real momenta we can use the identity, 



p+Vp^ = Vp+p- = \Ip-p+ = \lpi+pl , 



(9) 



to write the corresponding spinor. 



Aa(p) 



(p. 

[pi 



Pi + ip2 



Pl+pl 



■iff) 



(10) 

(11) 



2.2 Four- dimensional Spinors 



Positive and negative energy solutions of the four dimensional massless Dirac 
equation are identical up to normalization conventions. By closely following 
the definitions of [57], the solutions of definite hehcity 
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u±{k) = 1(1 ± -f,)u{k) and v^{k) = i(l ± 75)t;(fc) 
and their conjugates 



:i2) 



u±{k) = u{kmiTl5) and v^{k) = v{k)l{l T I5) ■ (13) 

can be chosen to be equal to each otheiQ- For the numerical evaluation of 
spinor products with slashed matrices insertion, we use the Dirac 7 matrices 
defined as, 
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1 

-1 



r 



/ 

V 



a' 
-a' 



75 



1 

1 



(14) 



where the entries and 1 are 2x2- matrices. In this representation, the massless 
spinors can be chosen as follows. 



u+{k) = v.{k) = -j= 



Up) 



u^{k) 



Up) 
-Up) 



(15) 



2.2.1 Label Representation 

The spinor-helicity formalism is based on the algebraic manipulation of basic 
objects such as the spinor inner products. They are independent of the ex- 
plicit representation used for expressing the spinors themselves. It is useful 
to have a notation for the spinors that is not bound to a particular explicit 
representation. In S@M we will call this notation the "Label" representation. 
We define 



u+{ki) = v^{ki) = \kf) =\i) = X{ki) , 

u.{h) = v+{ki) = \ki) = \i\ = \{ki), (16) 
and for the conjugate spinors 



u+{ki) = v^{ki) = 



kt 



u-{ki) = v+{ki) = (A;. I = {i\ = X{ki) 



(17) 



^ Note that for negative energy solutions, the helicity is the negative of the chiraUty 
or 75 eigenvalue 
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2.3 Spinor Products 



The scalar (or inner) products for massless spinors are defined as, 



{ij) = {kr\kp = u.{h)u+{k,) = \\h)\a{k,) = X{h)X{kj), (18) 
[ij] = {kt\kj) = ^I4^«_(%) = \a{h)\\ki) = X{ki)Ukj). (19) 

Tfie fielicity projection implies that products like [i\j) vanish. 

In the rest of the paper we will call spinor products the ( )- and []-type spinor 
inner products, and not the external spinor-product for the spinorial decom- 
position of slashed-matrices, see Eq.(7). 

The spinor products arc, up to a phase, square roots of Lorentz products, in 
fact they are related to momenta inner product through the identity, 

{ij)\Ji] = (r = IV(i(l - 75) i^j) = 2ki ■ kj = Sij, (20) 

where Sij = {ki + kj)"^ = 2ki ■ kj. 
We also have the useful identities: 
Gordon identity: 

=[i|7'^|i) = 2A;f. (21) 

Projection operator: 

|.^[.|^(l+75) |.]^.|^^75 + (22) 

Antisymmetry : 

{ji) = -{ij), [ji\ = -[iji {ii) = [H] = . (23) 
Schouten identity: 

{i3){kl) = {ik){jl) + {il){kj). (24) 
Spinor re-definition: 
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m-\m, m^-{k(i)\. (25) 



2.4 Massive Spinors Products 



In the spinor-hclicity formalism, spinors associated to particles of momentum 
P/which are solutions of the massive Dirac equation, 



u 



Hpi) {fj - mi) u^pi) = , 



can be as well represented with bra-ket notation: 



u+{pj)^\I), u+jpj) ^ {I\ , (26) 
u-ipi) = \I], u-ipj) = [I\, (27) 

where the angle-bracket denote the two spin-polarization states with respect 
to a fixed reference axis, rather than helicity states. Spinor products in this 
case can be written as [70, 71] 



u ipi)u^ipj) = {IJ) , u+{pj)u (pj) = [IJ] , 

^^^H^u^ipj) = [IJ) , ^F(^«-(pj) = [IJ] . (28) 

One can use the light-cone decomposition of a massive spinor in terms of two 
massless one [86]. Accordingly, by introducing an arbitrary massless reference 
momentum g, we can construct a massless momentum (pj) associated to any 
massive one (p/). 



P'-P'-^f-P'-^.^l'^ pNO = g^ (29) 
Zpi ■ q ■Zpi- q 

Correspondingly, the spinor variables read. 



= + (30) 

= N + (31) 

(32) 

With the above formulas we can express the spinor product for massive par- 
ticles of Eqs (28) in terms of the ones involving only massless spinors [71] 
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(/J) 
[IJ] 




(33) 
(34) 



[I J) 



(35) 



(7 J] 



(36) 



where Siq — {pi+qY — 2pi-q — [iq]{qi). We notice that in the r.h.s of the above 
equations only massless spinors do appear. Therefore the above definitions of 
spinor products associated to massive particles can have their straightforward 
implementation in SOM. 



3 Mathematica Implementation 

In the following section we present the Mathematica implementation SOM. The 
Mathematica notebook file S@M_Def initions .nb containing all the examples 
shown in this section is distributed with the package. 

3.1 Input and Output 

In the notebook, both the input and the output has been designed to look like 
the familiar bra-ket representation. 



In the console version, the functions can only by inputed using their names, 
in the notebook, several other input methods are available. 

• Output from a previous evaluation can be copy-pasted as input for the next 
evaluation. 

• Some basic functions and objects like the spinor products are available from 
a palette that is opened when the package is loaded within the notebook. 

• Spinor products can be entered using the characters sequnce "ESC < ESC" , 
"ESC > ESC","[", "]", "— " if the input format is TraditionalForm. 

The examples in the following sections are displayed as they appear in the 
Mathematica notebook. 

Similar output can be obtained in the console version by using the command 
TraditionalForm. The package is loaded using 



(•,•),[•,•]>(•>•>•]■■■ 
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<< Splnors ^ 

SPINORS @ MATHEMATICA (S@M) 

Version: S@M 1.0 (7-OCT-2007) 

Authors: 

Daniel Maitre (SLAC), 

Pierpaolo Mastrolia (University of Zurich) 

A list of all functions provided by the package 
is stored in the variable 
$Spinor Functions 



The package should be loaded at the beginning of the Mathematica session. 



3.2 General Structure 



Massless spinor variables are the fundamental objects for the construction of 
both inner spinor- products (associated to Lorentz invariants), and external 
spinor-products (associated to slashed-matrices). 

In SOM, massless spinors are expected to be labeled by either symbols or in- 
tegers. The spinor products are independent of the explicit representation of 
the spinors. Objects that do not refer to an exphcit representation, hke spinor 
products in their (•,•), [•,•] form, are said to be in the "label" representa- 
tion. To allow more flexibility in the use of the package two explicit spinor 
representations have been implemented, the two-dimensional (Weyl) and the 
four- dimensional (Dirac) ones. The different representations and the names 
are summarized in the following table. 



Label 


2 dimensional 


4 dimensional 






u+{ks) , u^{ks) 


s 


La[s],CLa[s] 


USpa[s], UbarSpaEs] 




\{ks),Uks) 


u-{ks), u+{ks) 


s 


Lat [s] ,CLat [s] 


USpb [s] , UbarSpb [s] 



The functions provided by the packages act on (defined) spinor-labels, and 
work in all the representations. The flexibility of S@M relies in the multiple 
use of a given symbol defined as massless spinor, which can represent at the 
same time either the spinor itself or its associated slashed matrix, with the 
automatic understanding of its interpretation according to the context. 
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3.3 New Functions 



The new functions introduced by the package are described in the following 
section. A list of all functions can be found in Appendix A. 

3.3.1 Spinors Declaration 

In S@M objects called (and declared as) spinors are considered to be the solution 
of the massless Dirac equation. That is not a restriction on the usability of the 
package, since solutions of the massive Dirac equation can be constructed from 
massless spinors. The symbol representing a given spinor is used to represent 
at the same time, and depending on the context, the spinor itself, its vector 
or the corresponding slashed matrix. 

DeclareSpinor 

The function DeclareSpinor can be called with one or a sequence of argu- 
ments. It declares its arguments to be spinors. If undeclared variables are 
used as spinors, some automatic properties will not be applied and most 
functions cannot be used. 



DeclareSpinor [a, b, c. 


d, e, f , i. 


j, k, 1, s, t, u] 


{a, b, c, d, e, f, i, j. 


k, 1, s, t. 


u} added to the list of spinors 



Integer labels for spinors do not have to be declared, for more details, see 
the section on Sp below. If a symbol is defined as a spinor, it can also be 
used to represent both its Lorentz vector (see page 14) or its corresponding 
slashed matrix. 



SpinorQ 

SpinorQ tests whether its argument has been declared as a spinor or not, 
it returns True if so and False otherwise. It can be used for example in 
patterns. 

SpinorQ [a] 

MatchQ[b, _? SpinorQ] 

True 

True 



UndeclareSpinor 

The function UndeclareSpinor removes its argument from the list of spinors. 

Sp 

Spinors can be labeled by integers using the function Sp. The object Sp [n] 
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is considered as a spinor. In the scalar products Spaa, Spab, Spba and Spbb 
described below, integer arguments are automatically wrapped into the Sp 
function. 

Sp[l] 

Sp[l] // FullForm 
SpinorQ[Sp[l] ] 

1 

Sp[l] 
True 

In StandardForm and TraditionalForm, the function Sp is not displayed, 
only its argument. 

The function can be made visible by using FullForm. 

3.3.2 Spinor Representations 

La, Lat, CLa, CLat 

The two-dimensional representations of the spinor s are given in the follow- 
ing table. 



A(A;,) 




Uks) 


Uks) 


La[s] 


Lat [s] 


CLaEs] 


CLat [s] 



The arguments of these functions are spinor labels or integers. In the latter 
case, the argument, say i, is automatically converted to the corresponding 
spinor label, using Sp[i] (see above). 

La[s] , Lat [s] , CLa[s] and CLat [s] are linear in their spinor argument. 

La[2 a - b] 

% / / TraditionalForm 

CLat [a + b + c] 

% // TraditionalForm 

2 La[a] -La[b] 
2X{a)-X(b) 

CLat [a] +CLat[b] +CLat[c] 
A(a) + A® + A(c) 

The contraction of two different two-dimensional spinors is implemented 
using the Mathematica Dot operator and automatically displayed in a stan- 
dard order (appropriate for the numerical evaluation, see page 21). 

CLa[l] .La[2] // TraditionalForm 
La [2]. CLa [1] //TraditionalForm 

A(1).A(2) 

A(l)..\(2) 



USpa, USpb, UbarSpa, UbarSpb 
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The four-dimensional representation of the spinor s are given in the follow- 
ing table. 



u+{ks) 


u-{ks) 


u-{ks) 


u+{ks) 


USpa[s] 


USpb[s] 


UbarSpa [s] 


UbarSpb[s] 



These four functions are linear in their spinor argument. 

USpa [2 a - b] 

% / / TraditionalForm 

UbarSpb [a + b + c] 

% // TraditionalForm 

2 USpa [a] - USpa[b] 
2M+(a) — M+(6) 

UbarSpb [a] + UbarSpb [b] + UbarSpb [c] 
+ Jmb) + «7(c) 

The contraction of two different two-dimensional spinors is implemented 
using the Mathematica Dot operator and automatically displayed in a stan- 
dard order (appropriate for the numerical evaluation, see page 21). 

UbarSpa [1] .USpa [2] //TraditionalForm 
USpa [2] .UbarSpa [1] //TraditionalForm 

Mr(i).H+(2) 
Mr(i).H+(2) 



3.3.3 Lorentz Vectors 

In SOM we call Lorentz vector any 4-dim vector of the form, k'^ = {k^, k^, k'^, k^). 



DeclareLVector 

The function DeclareLVector can be called with one or a sequence of ar- 
guments. It declares its arguments as Lorentz vectors. 

Momenta associated to spinors (declared through DeclareSpinor) do not 
need to be declared, and one can use the symbol of the spinor to represent 
the corresponding Lorentz vector as well. To represent the Lorentz vector 
of a spinor labeled by an integer i, one may use Sp [i] . 

DeclareLVector [pi , p2, p3] 

{pi, p2, p3} added to the list of Lorentz vectors 



UndeclareLVector 

The function UnDeclareLVector removes its argument from the list of the 
Lorentz vectors. 
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UndeclareLVector [p3 ] 

p3 removed from the list of Lorentz vectors 



LVectorQ 

LVectorQ tests whether its argument can be interpreted as a Lorentz vector 
or not. It returns True if so and False otherwise. It can be used for example 
in patterns. 

LVectorQ [pi] 
LVectorQ [p3] 
LVectorQ [Sp[ 5] ] 

True 

False 
True 

MatchQ[pl, _? LVectorQ] 

True 



3.3.4 Minkowski Products 
MP [p,q] 

The function MP [p , q] represents the Minkowski product 

00 - -> 
p-q^p q -p-q. 

The two arguments can be either spinors or Lorentz vectors or hnear com- 
binations thereof. 

Integer arguments are interpreted as the four vectors associated to the 
spinor Sp [i] . MP is symmetric in its arguments, so they are automatically 
sorted. The vector product is linear. 

MP[2, 1] 

MP[2pl, 3p2-5Sp[l]] 

MP[1, 2] 

6MP[pl, p2] -10MP[pl, 1] 



MP2 [p] 

The function MP2 [p] is a shortcut for MP [p,p] . 

MP2[3pl -p2] 

9MP[pl, pi] -6MP[pl, p2] +MP[p2, p2] 

The function MP can be used with explicit four vector representations: 
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MP[{5.3, 1.4, -1.2, 3.7), {1.2, -0.7, 0, 1.1}] 

3.27 



3.3.5 Invariants 



s[i, j] 

The function s [i , j] represents the kinematic invariant given by the square 
of the sum of two momenta, 

Sij = {Pi+Pjf- 

The two arguments can be either spinors or Lorentz vectors. 

Integer arguments are interpreted as the four vectors associated to the 
spinor Sp [i] . Since the scalar product s is symmetric in its arguments, 
they are automatically sorted. 

s[a, b] 
s[4, 1] 

Sab 
Slx4 

The function s also accepts more than two arguments (which can be spinors 
or Lorentz vectors) for multi-particle invariants, 

Si...j = {Pi + ... 



s[l, 2, 3] 
s[4, 3, 2] 

Slx2x3 
32x3x4 



3.3.6 Slashed Matrices 

Slashed matrices are in general contractions of Lorentz momenta with gamma- 
matrices ^ — P^7^. There are three representations for the slashed matrices 
in the package, as summarized in the following table. 
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Label 


2 dimensional 


4 dimensional 


f 


pda p 
^ •) ^ ail 




Sm[P] 


Sm2 [P] , CSm2 [P] 


Sm4 [P] 


m + \a)[h\ 




u^{a).U-{b) + U-{b).u^{a) 


SmBA [b,a] 


SmBA2 [b , a] , CSmBA2 [b , a] 


SmBA4 [b,a] 



In S@M slashed matrices are used either inside spinor products or as explicit 
matrices. There are four different types of slashed matrices. 



• slashed matrices corresponding to a declared (massless) spinor. These slashed 
matrices are labeled with the same symbol as their spinor. The slashed ma- 
trix associated with the massless spinor a are represented by: Sm [a] , in the 
label representation; Sm2 [a] or CSm2 [a] , in the two-dimensional represen- 
tation; and Sm4 [a] , in the four-dimensional representation. Inside spinor 
products in the label representation, the function Sm is not necessary and 
can be omitted. 

• external products of spinors. These slashed matrices are represented by the 
functions: SmBA, in the label representation; SmBA2 or CSmBA2, in the two 
dimensional representation; and SmBA4, in the four-dimensional representa- 
tion. 

• slashed matrices corresponding to a declared (possibly massive) vector. 
These slashed matrices are labeled with the same symbol as their vector. 

The slashed matrix associated with the vector P are represented by: Sm [P] in 
the label representation; Sm2 [P] or CSm2 [P] , in the two dimensional repre- 
sentation; and Sm4 [p] , in the four-dimensional representation. Inside spinor 
products in the label representation, the function Sm is not necessary and 
can be omitted. 

• other slashed matrices unrelated to a declared vector or spinor. These slashed 
matrices have to be declared. The slashed matrix declared with the symbol 
Q is represented by Sm [Q] in the label representation, Sm2 [Q] or CSm2 [Q] , 
in the two dimensional representation and Sm4 [Q] , in the four-dimensional 
representation. Inside spinor products in the label representation, the func- 
tion Sm is not necessary and can be omitted. 



DeclareSMatrix 

The function DeclareSMatrix can be called with one or a sequence of ar- 
guments. It declares its arguments to be slashed matrices. Slashed matrices 
corresponding to declared spinors and Lorentz vectors are arc labeled by 
the same symbol as the spinor or Lorentz vector and thus do not need to be 
re-defined. If undeclared variables are used as slashed matrices, some auto- 
matic properties will not be applied, and most functions cannot be used. 
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DeclareSMatrix[P, Q, R, S, T] 

{P, Q, R, S, T} added to the list of slashed matrices 



UndeclareSMatrix 

The function UndeclareSMatrix removes its argument from the hst of the 
slashed matrices. 

SMatrixQ 

SMatrixQ tests whether its argument has been declared as a slashed matrix 
or not, it returns True if so and False otherwise. It can be used for example 
in patterns. 

SMatrixQ [F] 
MatchQ[Q, _? SMatrixQ] 

True 

True 



SmBA 

The object SmBA [b, a] represents slashed matrices formed by the tensor 
product of two spinors, like 

\b]{a\ + \a)[b\ . 

The arguments a and b are spinors labels. SmBA is linear in both arguments. 
If the two arguments are equal, SmBA [a, a] is automatically replaced by 
Sm [a] . 

SinBA[l, 2] 
SniBA[b + c, a] 
SniBA[a, a] 

SmBA[l, 2] 

SmBA[b, a] +SmBA[c, a] 
Sm[a] 



3.3.7 Slashed Matrices Representations 
Sm 

The object Sm is used for slashed matrices corresponding to previously de- 
clared spinors and Lorentz vectors. Sm can be called with one argument, 

being either a spinor label or a vector label. In particular, slashed matrices 
associated either to spinors (declared through DeclareSpinor), or to vec- 
tors (declared through DeclareLVector) are automatically declared. One 
can use the symbol of the spinor, say s, or the one of the vector, say P 
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to represent the corresponding slashed matrix by means of Sm [s] or Sm [P] 

respectively. 

The object Sm is linear in its argument. 

SMatrlxQ[Sni[a] ] 
Sni[a + 2 b] 

True 

Sm[a] + 2 Sm[b] 



Sm2, CSm2 

The two-dimensional representations of the slashed matrix P are given in 
the following table. 



pita 


p . 


Sm2 [P] 


CSm2 [P] 



The functions Sm2 and CSm2 are linear in their spinor argument. The con- 
traction with another slashed matrix or with a spinor is implemented using 
the Mathematica Dot operator. 



Sni2 [a - b] 




Sni2 [a + b] .CSm2 [c + 


d] 


Sm2 [a] - Sm2 [b] 




Sm2 [a] .CSm2 [c] + Sm2 


[a] .CSm2 [d] + Sm2 [b] .CSm2 [c] + Sm2 [b] .CSm2 [d] 



Sin4 

The four-dimensional representations of the slashed matrix P is given by 
Sm4[P]. This function is linear in its spinor argument. The contraction of 
two different four-dimensional matrices or with four-dimensional spinors is 
implemented using the Mathematica Dot operator. 

Sm4[a - b] 

Sm4 [a] .Sm4 [a + b] 

Sm4 [a] - Sm4 [b] 
Sm4 [a] .Sm4 [b] 



SmBA2 

In the two-dimensional representation, the slashed matrices build from two 
spinors 



A .A" and \a\ 

are represented by the two function SmBA2 [b , a] , CSmBA2 [b , a] respectively. 
These functions are linear in both spinor arguments. 
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SniBA2[b, a + b] 
CSniBA2[b, a + Sp[l]] 

Sm2 [b] + SmBA2 [b, a] 
CSmBA2[b, a] + CSmBA2 [b, 1] 

Their multiplication with other slashed matrices are implemented using the 
Mathematica Dot operator. SinBA2 and CSinBA2 are replaced by their tensor 
product representation when they are inserted in a spinor chain. 

Lat[l] .SniBA2[2, 3] .La[4] // TraditionalForm 

CLa[a] .CSinBA2[b, c] .Sni2[d] .La[e] // TraditionalForm 

-A(3).A(4)A(2).A(1) 
A(a).A(c) A(6).Sm2(rf).A(e) 



SmBA4 

The slashed matrix build from two spinors, 

\b]{a\ + \b]{a\ , 

is represented in the four-dimensional representation by the function SmBA4 [b , a] . 
This function is linear in both spinor arguments. 

SniBA4[b, a + b] 

Sm4 [b] + SmBA4 [b, a] 

Its multiplication with other slashed matrices is implemented using the 
Mathematica Dot operator. SmBA4 is replaced by its tensor product rep- 
resentation when inserted in a spinor chain. 

ObarSpafl] .SniBA4[2, 3] .USpb[4] //TraditionalForm 
UbarSpb[a] .SinBA4[b, c] .Sm4[d] .USpb[e] // TraditionalForm 

-wr(l).»+(3)M7(4).M_(2) 
-itjTW-M-Ca) Hr(c).Sm4(rf).M_(e) 



3.3.8 Spinor Products 



Spinor products are represented in S@M by four different objects: Spaa, Spab, 
Spba and Spbb, according to the following table. 



{a...b) 


{a...b] 


[a...b) 


[a...b] 


Spaa [a, . . . ,b] 


Spab [a, . . . ,b] 


Spba [a, . . . ,b] 


Spbb [a, . . . ,b] 



The left- and right-most arguments are spinors and the intermediate argu- 
ments arc slashed matrices or objects that can be interpreted as slashed ma- 
trices such as spinors or Lorentz vectors. Spaa and Spbb expect an even num- 
ber of (declared) arguments whereas Spab and Spba expect an odd number 
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of (declared) arguments. If not all the arguments are integers (and therefore 
automatically interpreted as spinors) or declared either as spinors, Lorentz 
vectors or slashed matrices, the properties listed below cannot be applied. 

-Standeird order 

The spinor products have a normal ordering for their arguments. If the 
rightmost and leftmost elements are spinors, the middle elements are slashed 
matrices (or can be interpreted as such) and in addition if the spinors are 
not in the standard order, the spinor product is ordered using the identities 

(6a) = — {ah) , [6a] = — [a6] , 
{h\Q...P\a) = - {a\P...Q\h) , [b\Q...P\a] = - [a\P...Q\h] , 
[h\P\a) = {a\P\b] , [b\Q...P\a) = {a\P...Q\b] . 

A special case of these identities is the on-shell condition 

(a a) — 0, [a a] = . 



Spbb[a, a] 

Spaa [nospinor, nospinor] 
Spaa[b, a] 
Spbb[a, P, Q, b] 
Spaa [a, b] Spbb[a, b] 



(nospinor | nospinor} 

-{a I b) 

- lb \ Q \ P \ a] 
-<a \b) [b\ a] 

The normal ordering of the Spaa and Spbb products are opposite so that 
the products (a b) [b a] are displayed in this usual way. 

The same rules apply to the spinor products written in the two- and four 
dimensional representations. 

CLa[2] .La[l] 

Lat[3] .CLat[4] 

CLa[l] .CSin2[3] .Sin2[2] .La[l] 

-CLa[l] .La[2] 

-Lat[4] .CLat[3] 

-CLa[l] .CSm2[2] .Sm2[3] .La[l] 

CLat[3] .CSm2[2] .CLa[l] 

Lat[3] .Sm2[2] .La[l] 

UbarSpa[2] .USpa[l] 

UbarSpb[3] .USpb[4] 

UbarSpa[l] .Sm4[3] .Sm4[2] .USpa[l] 

-UbarSpa[l] .USpa[2] 

-UbarSpb[4] .USpb[3] 

-UbarSpa[l] .Sm4[2] .Sm4[3] .USpa[l] 

UbarSpb[3] .Sm4[2] .USpa[l] 

UbarSpa[l] .Sm4[2] .USpb[3] 
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-Linearity 

If the arguments of a spinor product have been defined as spinors using 
DeclareSpinor or DeclareSMatrix, then linear combinations of spinors 
are automatically expanded. 

Spaa[a + b, c] 

Spaa [2 b, c] 

Spab[a + 2b, P, 3 c - d] 

<a I c> + <Jb I c> 

2 <jb I c> 

3 {a \ P \ c] - (a \ P \ d] + 6 (b \ P \ c] - 2 {b \ P \ d] 
Lat [a + b] . Sin2 [2 c + d] . La [e] // TraditionalForm 
UbarSpb[a + b] .Sm4[2c + d] .USpa[e] // TraditionalForm 

2A(a).Sm2(c).A(e) + l(a).Sm2(d).X(e) + 2 A(6).Sm2(c).A(e) + kb).Sm2(d).Xie) 
2uI(e).Sm4(c).u-{a) + 2uI(e).Sm4(c).u-{b) + uI{e).Sm4(d).u-(,a) + Mr(e).Sm4(rf)-«-(6) 

The linearity works for integer labels too, but one has to write the Sp func- 
tion explicitly, since the sum or product of the integers is done before wrap- 
ping the result with the function Sp. 

Spaa[l + 2,2 coeff] 
Spaa[coeff Sp[l] , 2] 
Spaa[Sp[l] + Sp[2] , 3] 

<3 I 2 coeff) 
coeff <1 I 2> 
<1 I 3> + <2 I 3) 



-Synt£LK correction 

When all the arguments of a spinor product are declared either as spinor, 

Lorcntz vector or slashed matrix and their number docs not match the ex- 
pected number for the particular spinor product type (even for Spaa, Spbb, 
for odd for Spab, Spba), the type of the spinor product is changed automat- 
ically and a warning is issued. 



Spaa[l, Q, 2] 






















Spab[l, 2] 






















Spba[l, Q, P, 2] 






















Spbb[l, Q, 2] 






















Spaa: :wrongNbrSM 






















Wrong number of 


slashed 
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. > product 
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.] 
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Spab: :wrongNbrSM 






















Wrong number of 


slashed 


matrices 


in 


the <. 


. ] product 


Automatically 


changed 


to 


< . 


. > 


<1|2> 






















Spba: :wrongNbrSM 
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- [2 1 P 1 Q 1 1] 






















Spbb: :wrongNbrSM 
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-Slashed matrices insertion 
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The Dirac equation and on-shell condition 

ji \a) = 0, ^ |a] = 0, = Ci' = 0, 

are used for spinors used as slashed matrices. The inserted SmBA objects are 
automatically expanded, as shown in the following examples. 

Spab[l, a, a] 
Spbb[l, a, a, 2] 




Lat[a] .Sni2[a] .La[b] 

CLa[l] .CSm2[a] .Sni2[a] .La[2] 




UbarSpb[a] .Sni4[a] .USpa[b] 
UbarSpa[l] .Sm4[a] .Sni4[a] .USpa[2] 




Spab[l, SinBA[2, 3] , 4] 

Spaa[l, SinBA[2, 3], SniBA[4, 5], 6] 

-<1 I 3> [4 I 2] 

-<1 I 3) <5 I 6> [4 I 2] 

Lat[l] .SniBA2[b, a] .La[2] // TraditionalForm 
UbarSpa[l] .SniEA4 [b, a] .USpb[2] //TraditionalForm 

A(fl).A(2)A(l).A(fo) 

Mr(a).M+(l) im2).u-(b) 



3.3.9 Spinor Manipulations 

ExpeindSToSpinors , ConvertSpinorsToS 

The function ExpandSToSpinors, ConvertSpinorsToS convert invariants s 
to spinor products and conversely. 

ExpandSToSpinors [s[l, 2] s[2, 3]] 

<1 I 2> <2 I 3> [2 I 1] [3 I 2] 
ConvertSpinorsToS [%] 

Slx2 S2x3 

Lat[l] .Sin2[2] .CSni2[l] .Sin2[3] .La[l] 
ConvertSpinorsToS [%] 

Lat[l] .Sm2[2] .CSm2[l] .Sm2[3] .La[l] 

Slx2 Slx3 

UbarSpa[l] .Sni4[2] .Sin4[l] .USpa[3] 
ConvertSpinorsToS [%] // TraditionalForm 

UbarSpa[l] .Sm4[2] .Sm4[l] .USpa[3] 

Mr(l).H+(3) ilx2 



SpOpen, SpClose 

The function SpOpen decomposes spinor chains containing any slashed ma- 
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trix that corresponds to a massless spinor with products of smaller spinor 
chains, by applying the definition of such a matrix in terms of its opposite- 
chirality spinors, 

^ ^ \k]{k\ + \k)[k\ . 



The function SpClose has the reverse efi'ect as that of SpOpen. It attempts 
to replace products of spinor products with spinor chains containing slashed 
matrices. 

Both the functions can take either one or two arguments. The first argument 
is the expression to be manipulated; the second argument must be a spinor. 
With one argument, the functions open or close wherever possible. 

SpOpen[Spaa[l, 2, 3, 4, 5, 6]] 

<1 I 2> <3 I 4> <5 I 6> [3 I 2] [5 | 4] 

Spqpen[Spab[l, 2, P, 3, 4]] 

-<1 I 2> <3 I P I 2] [4 I 3] 

SpOpen [Spaa [1, 2, F, 3, 4, Q, 5, 6]] 

-<1 I 2> <5 I 6> <3 I P I 2] <4 I Q I 5] [4 I 3] 
SpPpen[Spaa[l, 2, 3, 4, 1, 2, 3, 4, 2, 3, 1, 2]] 
SpClose [%] 

<1 I 2)^ (2 I 3> <3 I 4>^ [3 I 1] [3 I 2]^ [4 | 1] [4 | 2] 
<1 12|3|4| 1 |2|4|3|2| 1 |3|2> 

Lat[5] .Sm2[4] .CSm2[3] .Sm2[2] .La[l] // TraditionalForm 
SpOpen [%] //TraditionalForm 
SpClose [%] //TraditionalForm 

A(5).Sm2(4).CSm2(3).Sm2(2).A(l ) 
A(1).A(2) A(3).A(4) A(3).A(2) A(5).A(4) 

A(5).Sm2(4).CSm2(3).Sm2(2).A(l) 

UbarSpa[l] .Sm4[2] .Sm4[3] .Sm4[4] .USpb[5] // TraditionalForm 
SpOpen [%] //TraditionalForm 
SpClose [%] //TraditionalForm 

Mr(l).Sm4(2).Sm4(3).Sm4(4).M_(5) 

«r(i).H+(2) Hr(3).M+(4) h7(3).m_(2) ii;:(5).M_(4) 
Mr(l).Sm4(2).Sm4(3).Sm4(4).M_(5) 

If there are different possibilities of reconstructing the spinor chain, SpClose 
does not search for the longest possible spinor chain. The result will depend 
on the ordering of the spinor labels and might not be invariant under rela- 
beling of the spinor labels. 

SpOpen [Spaa [1, 2, P, 4, Q, 2]] 
SpClose [%] 

<1 I 2) <2 I Q I 4] (4 I P I 2] 
<1 I 2> <4 I P I 2 I Q I 4] 

If a spinor is given as a second argument, SpOpen and SpClose will only 
open or close spinor chains containing this specified spinor. 
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SpPpen[Spat)[a, b, c, d, e] , c] 

(a \ b \ c] (c \ d \ e] 
SpOpen[Spab[a, b, c, d, e] ] 
SpClose[%, c] 

(a \ b) (c \ d> [c \ b] [e \ d] 

-(a \ b) (d \ c \ b] [e I d] 
Spaa[l, 2, 3, 4] // To2DiniSpinor 
Spppen[%, 2] 
SpClose [%] 

CLa[l] .CSm2[2] .Sm2[3] .La[4] 
CLa[l] .La[2] Lat[2] .Sm2[3] .La[4] 
CLa[l] .CSm2[2] .Sm2[3] .La[4] 

SpOpen[UbarSpb[d] .Sm4[c] .Sm4[b] .USpb[a]] // TraditionalForm 
SpClose [%, c] //TraditionalForm 

-Hr(6).H+(c) H7(6).M-(a) M7(rf)-M-(c) 

M7(fe).H-(a) Hr(6).Sm4(c).M_(rf) 



To2DimSpinor 

The function To2DimSpinor converts spinor products in the label represen- 
tation into the two-dimensional representation. 



To2DiinSpinor [Spaa[a, b, c, d] +Spab[a, c, d] ] ; 
% // TraditionalForm 

A(d).Sm2(c).A(a) + A(a).CSm2(6).Sm2(c).A(d) 
To2DimSpinor[Spab[a, F, d] ] ; 
% / / TraditionalForm 

A(rf).Sm2(P).A(a) 

There is an ambiguity in the conversion of explicit slashed matrices (see Sm), 
when they arc not embedded in a spinor chain. In those cases To2DimSpinor 
convert the slashed matrices to the Sm2 type. In case of Dot products of 
slashed matrices, the product is transformed to a Dot product starting with 

a. Sni2 matrix. 

To2DimSpinor[Sm[a + b] +SmBA[c, d] ] 

Sm2[a] +Sm2[b] +SmBA2[c, d] 
To2DimSpinor [Sm[a] .Sm[b] .Sm[c]] 

Sm2 [a] .CSm2 [b] . Sm2 [c] 



To4DimSpinor 

The function To4DimSpinor converts spinor products in the label represen- 
tation into the four-dimensional representation. 



To4DimSpinor[Spaa[a, b. 


c, d] +Spab[a, c, d] ] ; 


% // TraditionalForm 




M_(a).Sni4(c).M_(rf) + M_(a).Sm4(6).Sm4(c).H+(d) 


To4DimSpinor[Spab[a, F, 


d]]; 


% // TraditionalForm 




Mr(a).Sm4(/').M_(d) 




To4DimSpinor [Sm[a + b] + 


SmBA[c, d]] 


Sm4 [a] + Sm4 [b] + SmBA4 [c 


, d] 
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ToSpinorLabel 

The function To2SpinorLabel converts spinor products in the two- or four- 
dimensional representation into the label notation. 

Io2DiinSpinor [Spaa[a, b, c, d] +Spab[a, c, d] ] // TraditionalForm 
ToSpinorLabel [%] 

A(rf).Sm2(c).A(a) + A(a).CSm2(fe).Sm2(c).A(d) 

(a \ b \ c \ d) + (a \ c \ d] 

La[l].Lat[l] + La [2] .Lat [2] + La [1] . Lat [2] //TraditionalForm 
% / / ToSpinorLabel 

A(1).A{1) + A(1).A(2) + A(2).A{2) 

Sm[l] +Sm[2] +SmBA[2, 1] 

Sni2[F] .CSni2[Q] // ToSpinorLabel 

CSni2[P] .Sni2[Q] .La[3] // ToSpinorLabel 

P.Q 

(P.Q) [3] 

UbarSpb [ 4 ] . Sm4 [ 3 ] . Sm4 [ 2 ] . USpb [ 1 ] / / TraditionalForm 
% // ToSpinorLabel 

UbarSpa[2] .Sm4[R] .Sm4[Q] .Sm4[P] .USpb[l] // TraditionalForm 
% / / ToSpinorLabel 

M7(4).Sm4(3).Sm4(2).«_(l) 
[4|3|2|1] 

Mr(2).Sm4(R).Sm4(e).Sm4(/').H_(l) 

(2\R\Q\P\1] 



Compactif y 

Given a slashed matrix, say and a spinor of the |»)-type, say \b), one can 
construct the massless spinor 

which indeed is of the opposite chirality |«]-type. In fact, one can define, 

f>\b)^\P{b)], ^-[P(6)| , 

and, similarly, 

^\b]^\P{b)), [b\^^-{P{b)\ . 

In the package, the notation used for the spinor obtained by applying the 
slashed matrix to the spinor a is P [a] . 

SpinorQ[P[b]] 
SpinorQ[(P.Q) [b] ] 

True 

True 

These objects are recognized as spinors. The application of more than one 
slashed matrix to a spinor is done using the Mathematica Dot operator. The 
function Compactif y uses this definition to reduce spinor products with in- 
serted slashed matrices to spinor products of two spinor objects. 
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Compactify [Spab[a, P, b] ] 
Compactify [Spaa[a, P, Q, b] ] 
Compactify [Spbb[a, P, Q, b] ] 
Compactify [CLa [a] .CSin2[P] .Sin2[Q] .La[b]] 

wpm 

(a\(P.Qm) 
[(e.P)[a]|6] 

CLa[a] .La[ (P.Q) [b] ] 

One can specify a spinor as a second argument for Compactify. In this case 
the spinor products containing the given spinor are compactified in such a 
way that the specified spinor is left untouched. 

Compactify [Spab [a, P, b] +Spab[a, P, c] , b] 
Compactify [Spaa[a, P, Q, b] +Spab[a, P, c] , b] 

Compactify [UbarSpa[b] .Sm4[P] .USpb[c] +UbarSpa[a] .Sm4[P] .Sm4[Q] .USpa[b], b] // TraditionalForm 
Compactify [Lat[c] .Sm2[P] .La[b] +CLa[a] .CSm2[P] .Sm2[Q] .La[b], b] //TraditionalForm 

{a \ P \ c] - [P(a) I b] 

-{b I (Q.P) [a] > + <a I P I c] 

uZ(b).u^(P(c))--uZ(b).u4{Q.P)[a\) 

m.X(P{c))-mMQ.P)[a\) 

The functions ACompactif y [x, a] and BCompactif y [x, a] work the same 
way as Compactify [x, a] but onfy spinor products containing [a) or \a] re- 
spectively arc compactified. 

CoiiT'actify[Spab[l, P, 3] +Spab[3, P, 1], 3] 
ACompactif y [Spab [1, P, 3] +Spab[3, P, 1], 3] 
BCoii5>actify[Spab[l, P, 3] +Spab[3, P, 1], 3] 

-<P(1) I 3) + [3 I P(l) ] 
-<P(1) I 3> + <1 I P I 3] 
<3 I P I 1] + [3 I P(l)] 



UnCompact 

The function UnCompact uncompactifies the spinor products compactified 
with Compactify. 

Con^aactify [Spab[a, P, b] ] 
UnCompact [%] 

Compactify [Spab[a, P, Q, R, b] , b] 
UnCon5>act[%] 

Con5>actify[Lat[b] .Sm2[P] .La[a], b] 
UnCoii5>act [%] 

{a\Pm 
(a\P\b\ 

-[(R.Q.P) [a] I b] 
(a\P\Q\R\b] 
-Lat[P[a] ] .CLat[b] 
Lat[b] .Sm2[P] .La[a] 

One can specify a spinor as a second argument. In this case only the spinor 
products where the Dirac matrices are compactified onto the specified spinor 
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will be uncompactified. 



Conipactify[Spaa[a, P, Q, b] +Spab[b, F, c] , b] 
UnCompact [%] 
UnCompact[%%, c] 
UnCompact [%%%, a] 

(b I P{c) >-<£.! (Q.P) [a]) 
(a \ P \ Q \ b} + (b \ P \ c] 
-(b I (Q.P) [a]) + (b \ P \ c] 
(b \ P (c) ) + (a \ P \ Q \ b) 



Schouten 

The function Schouten applies the Schouten identities 

{ij){ke) - {ie){kj) + {ik){je) , [ij][ke] = [ie][kj] + [ik][je] . 

There are three different applications of the function. 
Schouten [x,i, j ,k,l] 

The function with four spinor arguments will search x for occurrences of 
the products {ij){k£) or [i and replace it using the above identities. 

Spaa[i, j] Spaa[k, 1] 
Schouten[%, i, j, k, 1] 

<i \ j) {k\ 1) 

-<i I J> <j I + <i I k) (j I J> 

CLa[i] .La[j] CLa[k] .La[l] // TraditionalForm 
Schouten[%, i, j, k, 1] // TraditionalForm 

X(i}.Xik) MJ}.Ml} - A(i).A(/) MJ).X(k) 



Schouten [x , i , j , k] 

The function with three spinor arguments will search for occurrences of 
the spinor product (ij) or [ij] and will try to use the Schouten identity 
to combine it with the spinor k. 



Spaa[i, j] Spaa[k, 1] Spaa[a, b] 


4'Spaa[i, j] Spaa[a, 1] Spaa[k, b] 


-<a 1 J> <£. 1 ;c> (1\ j) + (a\ b) (1 


j) (k 1 J> 


Schouten[%, i, j, k] 




-(a 1 J> ((b \ jy (1 \ k) - (b \ i> <j 


k)) +{a\b) (-<i 1 J> <j 1 i> + <i 1 k) (j 1 J>) 


Schouten [x, 1] 




The function with 


one spinor arguments will search for structures like 




{iu) [iu] 




{£s){£t) ' [es][et] ' 



and will use the Schouten identities to split them into partial fractions, 

{(u) _ {su) {tu) [iu] _ [su] [tu] 

{es){et) ~ {£s){st) ^ {£t){st) ' [is][it] ~ [is][st] ~ [it][st] ' 
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The function also works for spinor products with embedded Dirac matri- 
ces. 



Spaa[l, u] 




Spaa[l, s] Slpaa[l, 


t] 


Schouten[%, 1] 




(1 1 u) 




(1 \ s) (1\ t> 




<s 1 u> 


<t 1 u> 


<J \ s) (s\ t> <J 


t> <s 1 t> 



ASchouten, BSchouten 

The function ASchouten and BSchouten behave hkc Schouten but apply 
the Schouten identity selectively on |£)-variables and |£]-variables respec- 
tively. 



Spaa[l, u] 






Spbb[l, u] 










Spaa[l, S, T, s] Spaa[l, t] 


Spbb[l, s] Spbb[l, t] 










<i 1 u) 


[u 


1] 












(1 \ t) (1 \ S \ T \ s) 
Schouten [%, 1] 
ASchouten [%%, 1] 
BSchouten [%%%, 1] 

<t 1 u> 


[s 1 1] 


[t 
<s 


i] 

r 1 s 1 u> 


[u 


s] 




[u 1 t] 


(1 \ t) (s \ T \ S \ t) 
(t 1 u> 


<J 1 S 


T 
<s 


s> <s 1 r 1 s 1 t> 
r 1 s 1 u) 


[s 1 I] 
[u 


[t 

1] 


s] [t 


I] [t 1 s] 


<i 1 t) <s 1 r 1 s 1 t) 

<i 1 u) 


(1 1 S 
[u 


r 

s] 


s> <s 1 r 1 s 1 t> 

[u 1 t] 


[s 1 I] 


It 


i] 




<J 1 t> <J 1 s 1 r 1 s> 


[s 1 1] 


[t 


s] [t 1 J] [t 1 s] 











ASpinorReplace [x , a , n] , BSpinorReplace [x , a , n] 

The functions ASpinorReplace and BSpinorReplace replace the spinor a 
in expression x with n. ASpinorReplace only replaces the spinors \a) and 
BSpinorReplace replaces only \a]. Slashed matrices corresponding to the 
spinor a will be split according to 



^ = \a) [a\ + \a\ {a\ 
and the appropriate component will be replaced. 



ASpinorReplace [Spab [2, 1, 2], 2, a] 
BSpinorReplace [Spab [2, 1, 2], 2, a] 

<«|1|2] 
<2|l|a] 

CLa[a] .La[b] +Lat[a] .CLat[b] // TradltionalForm 
ASpinorReplace [%, b, c] //TradltionalForm 
DbarSpa[a] .USpa[b] -UbarSpb[b] .USpb[a] // TradltionalForm 
BSpinorReplace [%, b, c] // TradltionalForm 

M.a).X(b) - kb).ka) 

A(a).A(c) - A(6).A(a) 

uHa).u+(b) - imb).u-(a) 

uZ{a).u^(b) - imc).u-{a) 
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ASpinorShift [x,a,s] ,BSpinorShift [x,a,s] 

The functions ASpinorShift and BSpinorShift shift the spinor variable 
a in expression x with s. ASpinorShift only shifts the spinor \a) and 
BSpinorShift only shifts the spinors \a\. Slashed matrices corresponding 
to the spinor a will be split according to 

^ = \a) [a\ + \a] {a\ 

and the appropriate component will be shifted. 

ASpinorShift [Spaa [1, 2] Spbb[4, 2], 2, 3] 
BSpinorShift [Spaa [1, 2] Spbb[4, 2], 2, 3] 

(<1 I 2> + (1 I 3>) [4 I 2] 
<1 I 2> ([4 I 2] + [4 I 3]) 

ASpinorShift [CLa [a] .La[b] +Lat[a] .CLat[b], b, z c] // TraditionalForm 
BSpinorShift [UbarSpa [a] .USpa[b] -UbarSpb[b] .USpb[a] , b, z c] // TraditionalForm 

X(a).X{b) + zX{a).X(c) - A(6).A(a) 

iiZUi).ii+(h) -(77(/')./(-(") - z'ii^(c).u-(o) 

The shift parameter s will be interpreted as a spinor of the appropriate chi- 

rality. It can he a sum of s]:)inors. 

ASpinorShift [Spab [a, c, b] , c, xd+5e] 

(<a I c> + X <a I d> + 5 <a I e>) [c \ b] 

To account for more generic spinor definitions and shifts composite-spinors 
may be required. We can use 

a, P[b], (P.Q)[c], (P.Q.R)[d], ... 

as the shift argument of the Shift-functions to represent the compactified 
objects hke 

|a), |P(6)), |P.g(c)), \RQ.R{d)), ■■■ 
in the case of an A-shift. The expressions obtained using this kind of argu- 
ments can be later un-compactified by using UnCompact. 



DeclareSMatrix[M, Ml 


, M2] 




ASpinorShift [ Spaa [ a , 


b], b. 


M[c]] 


% / / UnCoii5>act 






{M, Ml, M2} added to 


the 11 


st of slashed matrices 


{a \ b) + (a \ M (c) ) 






(a \ b) + (a \ M \ c] 






ASpinorShift [ Spaa [ a , 


b], b. 


(M1.M2) [c] ] 


% / / UnCompact 






(a \ b) + (a \ (M1.M2) [ 


c]> 




<a 1 i) + <a 1 Ml 1 M2 1 c> 





ShiftBA[b,a,z] [x] 

The function ShiftBA[b,a,z] performs the shifts combination 
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\b]^ \b]- z\a\, \a) ^\a) + z\b). (37) 
The arguments a and b must be declared as spinors. 



ShiftBA[3, 2, z][Spaa[2, 4]] 
ShiftB&[3, 2, z][Spbb[3, 1]] 

<2 I 4> + z <3 I 4> 
-z [2 I 1] + [3 I 1] 

ShiftBA[c, b, z] [UbarSpa[a] .USpa[b]] // TraditionalForm 
ShiftBA[3, 2, z] [Lat[a] .Sm2[2] .La[b]] //TraditionalForm 

uZ(a).u^(b) + z Hr(a).M+(c) 

z A(6).A(3) A(2)i(a) + A(a).Sm2(2).A(6) 



3.3.10 Constants 

GammaO, Gammal, Gaimna2, Gamma3, GammaS are the 7-matrices in the rep- 
resentation (14). 

{GammaO, Gammal, Gamma2, Gamma3, Gamma5} // TraditionalForm 

SequenceFoon@@MatrixForm /@ {GammaO, Gammal, Gamma2, Gamma3, GammaS} //TraditionalForm 



{70,71.72,73,75} 



n 





















1 \ 


IQ 








-i ' 


10 





1 


\ 


10 





1 


\ 





1 
















1 











i 














-1 











1 








-1 










-1 











i 








-1 











1 











,0 








-1 




,-1 








, 










, 


.0 


1 





, 


,0 


1 





, 



The 7-matrices can be used as slashed matrices too. 



Spab[l, Gammal, 1] /2 

1 

- <1 I Ti I 1] 



ProjPlus, ProjMinus are the helicity projectors (1 ± 75)/2. 



{ProjPlus, ProjMinus} // TraditionalForm 

SequenceForm@@MatrixForm /@ {ProjPlus, ProjMinus} //TraditionalForm 


[P^, P-] 














1 

2 









i i 







1 

2 






i i 




1 

2 









lo i ij 







1 

2 ) 





The symbol $SpinorFunctions contains a list of all functions of the package. 



$SpinorsFunctions 

{ACompactif y, ASchouten, ASpinorReplace, ASpinorShif t, BCompactif y, BSchouten, BSpinorReplace, 
BSpinorShift, CLa, CLat, Compactify, ConvertSpinorsToS, CSm2, CSmBA2, DeclareLVector, 
DeclareLVectorMomentum, DeclareSMatrix, DeclareSpinor, DeclareSpinorMomentum, 
ExpandSToSpinors , GammaO, Gammal, Gamma2, GammaS, GammaS, GenMomenta, La, Lat, 
LVectorQ, MP, MP2, Num4V, Pf romBiSpinor, PfromCSm2, PfromSm2, PfromSm4, ProjMinus, 
ProjPlus, s, Schouten, ShiftBA, Sm, Sm2, Sm4, SMatrixQ, SmBA, SmBA2, SmBA4, Sp, Spaa, 
Spab, Spba, Spbb, SpClose, SpinorQ, SpOpen, To2DimSpinor, To4DimSpinor, UbarSpa, 
UbarSpb, UnCompact, UndeclareLVector, UndeclareSMatrix, UndeclareSpinor, USpa, USpb} 
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3.3.11 Numerics 



The spinor products have a numerical implementation. The first step to use 
this numerical implementation is to define or generate the four momenta for 
the spinors. For this there are several possibihties described in the following 
section. 



3.4 Momentum Generation 

DeclareSpinorMomentum 

This function takes as first argument the spinor whose momentum should 

be set, the second specifies the momentum. It can be one of the following. 

• The four vector in the form of a list {E,pl ,p2,p3}, (note that it is the 
responsibility of the user to make sure that the vector is indeed an on-shell 

vector) 

• two explicit spinors, the first in the A form and the second in the A form 
of (9). This can be used with the La and Lat functions. If the momentum 
associated with a spinor is defined using two spinors, these spinors will be 
used for the numerical evaluation of the spinor products. 

If the first argument was not declared as a spinor using DeclareSpinor, it 
will be declared. 



DeclareSpinorMomentTjm|^a, |V 14 , 1, 2, sjj 




DeclareSpinorMomentum l^b, |V 11 , 1, 1, sjj 




DeclareSpinorMomentum [c, {3, 2, 1, -2}] 

DeclareSpinorMbmentum[d, {{0.9}, {-0 . 5 - . 2 i} } , {{0.9, -0.5 + 
DeclareSpinorMomentum [e. La [a] //N, Lat[b] //N] 


2i}}] 


Momentum for spinor a set to |V 14 , 1, 2, sj. 




Momentum for spinor b set to |V 11 , 1, 1, 3|. 




Momentum for spinor c set to {3, 2, 1, -2}. 




Momentum for spinor d set to {0.55 + 0. i, -0.45 + 0. i, -0.18 + 


1, 0.26 + 0. 1} . 


Momentum for spinor e set to 
{3.4927 + 0.0766203 1, 1.00053 + 0.451416 1, 1.48451 + 0.0325661 1 


3.03298 - 0.0766203 1} . 



DeclareLVectorMomentum 

One can also define four- vectors that are not associated with a spinor (for ex- 
ample when they are not light-like) with the function DeclareLVectorMomentum. 
It takes two arguments, first the symbol for the vector to be defined and 
second the value of the four vector in the form of a fist. 



DeclareLVectorMbmentum[p, {1, 0, 0, 0}] 

{p} added to the list of Lorentz vectors 
Four Momentum p set to (1, 0, 0, 0}. 



GenM omenta 

The function GenMomenta[{sl , . . . ,sn}] generates arbitrary on-shell four 
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momenta for the spinors si, . . . ,sn. They are generated so that they sum 
up to zero momentum. 



GeiiManienta[{l, 2, 3, 4, d, e}] 

Momenta for the spinors d, e, 1, 2, 3, 4 generated. 



The function GenMoinenta[{sl, . . . ,sn}->{pO,pl ,p2,p3}] generates arbi- 
trary on-shell four momenta for the spinors si, . . . ,sn, so that the sum of 
these momenta is equal to the vector p — (po,Pi,P2,P3)- 



GeiiManienta[{l, 2, 3, 4} -> {2, 1, 


0, 


0)] 


Momenta for the spinors 1, 2, 


3, 


4 generated. 


GenMonienta|^ {a, b, c, d} -> j'/i". 


1, 


0, 1}] 


Momenta for the spinors a, b. 


c. 


d generated. 


GeiiMonienta[{l, 2, 3, 4} -> {0, 0, 


0, 


1)] 


Momenta for the spinors 1, 2, 


3, 


4 generated. 



The momenta can be produced with GenMomenta in a reproducible way 
by seeding the random number generator with the Mathematica command 
SeedRandom. 



SeedRandom [ 1234 ] 
GenMomenta [{1, 2, 3, 4, 5}] 

Momenta for the spinors 1, 2, 3, 4, 5 generated. 

GenMomenta can take as a second argument the precision with which the mo- 
menta should be generated. Too large precision slows down the numerical 
evaluation. 

GenMomenta [{1, 2, 3, 4, 5}, 50] 
Num4V[l] 

Momenta for the spinors 1, 2, 3, 4, 5 generated. 
{0.6081332974505249838643965958522341290362790390318372, 

-0.5208126293248651011171457745499236791066537010532711, 

-0.0332051403140836824289976445149074193363821989781881, 

0.3122142393621587 778454886924537509139639228243917401} 



3.5 Numerical Functions 



Once the momenta associated with the spinors are declared or generated, 
several numerical functions are accessible. 

Num4V 

The value of the momentum associated to the spinor a is accessible through 
the function Num4V. 

Num4V[a] 

{1.7018918588514296, 1.4060365615683739, -0.4260134319614734, 0.8590748759835581} 
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Num4V is linear in the same way as MP. 



Num4V[2a + b] 

{3.721558979299065, 2.908093303275366, -0.744906615859027, 2.001498467057566} 

The exphcit numerical representation of the spinors are accessible by ap- 
plying the Mathematica function N on the two- or four-dimensional repre- 
sentation, both are represented in a matrix notation. La, CLat, USpa and 
USpb are column vectors, whereas CLa, Lat, UbarSpa and UbarSpb are row 

vectors. 

{La[c] , CLa[c] } // TraditionalForm 
% / / N / / TraditionalForm 

{m,Mc)} 

I , ( 1 .38665 - 0.942001 i -2.48327 ) \ 

lU .38665- 0.942001 E j / 

{Lat [c] , CLat [c] } // TraditionalForm 

% / / N / / TraditionalForm 

{kc),Mc)} 

I H. 38665 + 0.942001 A, 

U 2.48327 1.38665 + 0.942001 i), | 

I 1-2.48327 )l 

USpa[c] // TraditionalForm 

% //N 

UbarSpa[c] // TraditionalForm 

% // N 

M+(C) 

{{1.75593}, {0.980508 - 0.666095 1}, {1.75593}, {0.980508 - 0.666095 1}} 
{{0.980508 - 0.666095 1, -1. 75593, 0.980508 - 0.666095 1, -1. 75593}} 

Spinor products with or without inserted slashed matrices corresponding to 
spinor momenta as well as the invariants s [a, . . . ,b] can be evaluated nu- 
merically. 

Spaa[a, b] // N 
Spbb[a, b] // N 
Spab[l, 3, 2] //N 
Spaa[l, 3, 4, 2] //N 
s[a, e] //N 

0.483013 - 0.427498 1 
-0.483013 - 0.427498 1 
1.43884 - 1.05022 i 
-1.64152 - 2.15347 1 
-0.909762 

Once the numerical value of four vectors arc defined, one can use them in 
spinor products with slashed matrices, vector products MP [p , q] , MP2 [p] 
and invariants s [p , q] . 
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DeclareLVector^mentum[p, 


{1, 


0, 


0, 0}] 


DeclareLVectorMbmentum[q, 


{1, 


0, 


0, -1}] 


Four Momentum p set to 


{1 


, 0, 


0, 


0}. 


{q} added to the list 


of 


^orent 


z vectors 


Four Momentum q set to 


{1 


, 0, 


0, 


1 \ 












Spaa [a, p, q, b] // N 




















1.36241 - 0.412794 i 










S[p, q] //N 










MP[p, Sp[2]] //N 










MP2[p + q-Sp[4]] //N 
3. 










0.725326 










8.76299 











PfrmSin2, PfroinCSm2 

The four vector can be extracted from the numerical two-dimensional repre- 
sentation of its slashed matrix using the function Pf roniSm2 or Pf romCSm2, 
depending on whether the matrix is of the Sm2 or CSm2 type. 



FfroniSin2[Sni2[2] //N] //Chop 




Num4V[2] // N 




PfromSin2[A(2) .A(l) //n] //Chop 




PfromCSin2[CSniBA2[2, 1] //N] //Chop 




{0.725326, -0.299499, -0.42836, -0. 


502897} 


{0.725326, -0.299499, -0.42836, -0. 


502897} 


{0.414319 - 0.235551 1, -0.432629 + 


.42751 1, -0.443834 - 0.176593 1, 0.0381312 + 0.235551 1} 


-'0.414319 - 0.235551 i, -0.432629 i 


.42751 1, -0.443834 - 0.176593 1, 0.0381312 i 0.235551 1} 


PfrmSm4 



The four vector can be extracted from the numerical four-dimensional rep- 
resentation of its slashed matrix using the function Pf romSm4. 



PfroniSin4[Sni4[3] //N] //Chop 




Num4V[3] // N 




{0.797155, 0.778221, 0.128146, 


115789} 


{0.797155, 0.778221, 0.128146, 


115789} 



4 Simple Examples 

In this section we will show three simple examples of the use of the S@M 

package. In the first example, we will re-derive the five gluon MHV amplitude 
^tree(^]^-^2-,3+,4+,5+) using the BCFW construction [67], in the second, 
we compute a box integral coefficient numerically using the quadruple cut 
technique. The third example illustrates the evaluation of cut integrals using 
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Fig. 1. BCFW construction of the five-gluon MHV amplitude using a |23) shift. 

the spinor integration method. The examples of this section are distributed 
with the package. 



4-1 BCFW Construction 



In this example, we will re-derive the very well-known five-gluon MHV ampli- 
tude ^^^^^(1", 2", 3+, 4+, 5+) using the BCFW construction [67]. The Mathe- 
matica notebook file S@M_BCFW . nb containing this example is distributed with 
the package. We start by defining the MHV amplitudes. 

MHV[a_, b_, c ] :=iSpaa[a, b] 4 / Times @@ Spaa @@@Partition[ {a, b, c, a}, 2, 1] 

MHVbar[a_, b_, c ] : = 

(-1) ^Length[{a, b, c}] iSpbb[a, b] ^ 4 / Times @@ Spbb @@@ Partition [ {a, b, c, a), 2, 1] 



Since we know the answer we can use it to check the calculation 



AO = MHV[1, 


2, 3, 4, 


5] 




i <1 1 2>" 




<1 1 5> <2 


3> <3 


4) <4 1 5> 



We will use the numerical implementation to check the result, so we need to 
generate a set of five on-shell momenta which sum up to zero momentum. 



GenMomenta[{l, 2, 3, 4, 5}] 

Momenta for the spinors 1, 2, 3, 4, 5 generated. 



We will use the shift 

A2 — * A2, A3 — s> A3 + zXs, A2 ^ A2 — ^A3, A3 A3. 



Now we are ready to start the computation. With the shift we chose, there is 
only one momentum partition contributing namely {5, 1,2}{3,4}, see Figure 
1. The amplitude is then given by 

DeclareSpinor [Q, mQ] 

{Q, mQ} added to the list of spinors 
1 

A = MHV[1, 2, Q, 5] MHVbar[3, 4, mQ] 

s34 

i <1 I 2)^ [4 I 3]^ 
s34 (Q \ 2) (Q\ 5> <1 I 5> [3 I mQ] [4 | mQ] 
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First we need to find the value of z that puts the propagator 



on-shell. 



zsol = Solve [ShiftBA[2, 3, z][s[3, 4]] ~ 0, z] // ExpandSToSpinors // SpOpen // Flatten 



<2 I 4> J 

We can check numerically that the formula for the amplitude is right, but for 
that we need a numerical expression for the spinors and \Q) associated 



with the shifted momentum Q = P3{z) + p4{z). This 


can be done with the 


functions DeclareSpinorMomentum and Pf ormSm2 




ShiftBA[2, 3, z] [Sin2[3] + Sm2[4] ] 




Qs = % / . zsol 




z CLat [3] .CLa[2] + Sm2 [3] + Sm2 [4] 




CLat[3] .CLa[2] <3 | 4> 

Sm2 [ 3 ] + Sm2 [ 4 ] 

(2 1 4) 




DeclareSpinorMbmentuni[Q, Ff roniSin2 [Qs // N] ] 




Momentum for spinor Q set to 




{-0.357954 - 0.401137 i, -0 . 032426 + . 0513217 i, 0.241304 - 0.486258 1 


, -0.588384 - 0.446289 1} . 


DeclareSpinorMomentum [mQ, -PfromSm2 [Qs // N] ] 




Momentum for spinor mQ set to 




{0.357954 + 0.401137 1, 0.032426 - 0.0513217 1, -0.241304 + 0.486258 1, 


0.588384 + 0.446289 1} . 



Now we can numerically evaluate the shifted amplitude and check that we get 
the right result. 



ShiftBA[2, 3, z] [A] 

i <1 I 2>^ [4 I 3]^ 
s34 (Q I 2> <Q I 5> <1 I 5> [3 I mQ] [4 | mQ] 
% /. s34 -» s[3, 4] // N 
-0.00573036 - 0.0112245 1 
N[AO] 

-0.00573036 - 0.0112245 1 

Here we see that the numerical result matches with the known MHV result. 
We can also recover the analytic formula for the MHV amplitude. For that we 
first need to convert the spinors \—Q) into \Q) using 

\-Q) = ±t\Q), \-Q] = ±i\Q] 



where the ± sign convention has no impact here since the spinor \Q] appears 
twice. 
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Fig. 2. Quadruple-cut for the determination of the coefficient of the box integral 
Ia{0, 0, 0, S45) in the one-loop five-gluon amplitude, with a gluon propagating in the 
loop. 



A / . mQ -> i Q 

i <1 I 2>^ [4 I 3]^ 
s34 (Q I 2> <Q I 5> <1 I 5> [3 I Q] [4 | Q] 
SpClose[%, Q] 

i <1 I 2>^ [4 I 3]^ 
s34 (1 I 5> <2 I Q I 3] <5 I Q I 4] 
ShiftBA[2, 3, z] [% /. Q^Sm[3] + Sm[4] ] /. zsol 
i <1 I 2>^ [4 I 3]^ 
s34 <1 I 5) (- ^""^il"'" + <3 I 5>) <2 I 4 I 3] 
% / . s34 -» s [3, 4] // ExpandSToSpinors // SpOpen // Sin5>lify 

i (1 I 2f 

<1 I 5> <3 I 4) (<2 I 5> <3 I 4> - <2 I 4> <3 | 5>) 

We can make the result look more familiar using the Schouten identity 

Schouten[%, 2, 5, 3, 4] 

i (1 I 2)' 

<1 I 5> <2 I 3> <3 I 4) <4 I 5> 

we can now check that this is the expected known result. 

MatchQ[%, MHV[1, 2, 3, 4, 5]] 

True 



4-2 A Box Coefficient 



We will compute the coefficient of the one mass box li{0, 0, 0, S45) of the one- 
loop five-gluon amplitude AI~^°°^{1~, 2~ , 3"*", 4"*', 5^) with a gluon propagating 
in the loop, see Figure 2. The Mathematica notebook file S(§M_Cut .nb contain- 
ing this example is distributed with the package. First we generate a momen- 
tum configuration with five momenta and define the loop momentum L. 
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SeedRandom [ 1 1 1 1 ] 
SenMoinenta[{l, 2, 3, 4, 5}] 

Momenta for the spinors 1, 2, 3, 4, 5 generated. 
DeclareLVectorMbnientum[L, {LO, LI, L2, L3}] 

{L} added to the list of Lorentz vectors 
Four Momentum L set to {LO, LI, L2, L3}. 

Note that using the command SeedRandom [ . . . ] before GenMomenta allows us 
to reproduce the same numerical momentum configuration in another session, 
so that comparison of numerical results are easier. 



Now we define the four propagator momenta 11,12,13,14 



11 = L 

12 = L - Sp[l] 

13 = 12 - Sp[2] 

14 = 13-Sp[3] 

L 

L - 1 

L - 1 - 2 

L - 1 - 2 - 3 



and solve for the components of the loop momenta L using the on-shell con- 
straints for the four propagators. 



sols = Solve [ { 






MP2[11] == 0, 






MP2[12] == 0, 






MP2[13] == 0, 






MF2[14] ~ 






} //N, {LO, LI, L2, L3}] 






{{LO ^ 1.12013 - 0.261024 i, LI ^ 


.492206 


- 1.6121 i. 


L2 ^ -1.84975 + 1. 77393 i, L3 -> - 


2.39887 


-1.576 75 i}, {LO^l. 12013 + 0.261024i. 


LI ^ 0.492206 + 1.6121 i, L2 ^ -1 


.84975 - 


1.77393 1, L3 ^ -2.39887 + 1.57675 i}) 



Using the first solution, we define the spinors corresponding to the propagator 
momenta (both incoming and outgoing). 



DeclareSpinorMomentuin [pll , Nuin4V[ll] /. sols[[l]]] 
DeclareSpinorMomentum[mll, -Num4V[ll] /. sols[[l]]] 
DeclareSpinorMomentuin [pl2 , Nuin4V[12] /. sols[[l]]] 
DeclareSpinorMomentuin [ml2, -Num4V[12] /. sols[[l]]] 
DeclareSpinorMomentum [pl3 , Num4V[13] /. sols[[l]]] 
DeclareSpinorMbmentum[ml3, -Num4V[13] /. sols[[l]]] 
DeclareSpinorMomentum [pl4 , Num4V[14] /. sols[[l]]] 
DeclareSpinorMbmentum[ml4, -Num4V[14] /. sols[[l]]] 




{pll} added to the list of spinors 
Momentum for spinor pll set to 
{1.12013 - 0.261024 i, 0.492206 - 1.6121 1, -1.84975 + 1.77393 1, -2.39S 

{mil} added to the list of spinors 


J87 - 1.57675 1} . 



The output has been shortened for readability. We can now compute the coef- 
ficient by inserting these spinors into the amplitudes at each corner of the box. 
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sollhl = MIIV[nill, pl4, 4, 5] MIIV[1, nil2, pll] MHVbar[nil3, pl2, 2] MHVbar[3, nil4, pl3] // N // Chop 



sollh2 = MHV[mll, pl4, 4, 5] MHVbar [nil2, pll, 1] MHV[pl2, 2, mlS] MHVbar[3, nil4, pl3] // N // Chop 

-41.0671 + 62.7569 i 



Where solhl and solh2 are the two contributions corresponding to the two 
(not obviously vanishing) hehcity configurations for the propagator momenta. 
The first hehcity configuration gives a vanishing coefficient. The second solu- 
tion for the loop momentum yields only vanishing solutions: 

DeclareSpinorMbnientum[pll, Nvmi4V[ll] /. sols [[2]]] 
DeclareSpinorMbmentum[nill, -Num4V[ll] /. sols[[2]]] 
DeclareSpinorMbnientum[pl2, Nvmi4V[12] /. sols [[2]]] 
DeclareSpinorMbmentum[nil2, -Num4V[12] /. sols[[2]]] 
DeclareSpinorMbnientum[pl3, Nvmi4V[13] /. sols [[2]]] 
DeclareSpinorMbmentum[nil3, -Num4V[13] /. sols[[2]]] 
DeclareSpinorMbnientum[pl4, Nvmi4V[14] /. sols [[2]]] 
DeclareSpinorMbmenti:im[ml4, -Num4V[14] /. sols[[2]]] 

Momentum for spinor pll set to 

{1.12013 + 0.261024 1, 0.492206 + 1.6121 1, -1.84975 - 1.77393 1, -2.39887 + 1.57675 1}. 
Momentum for spinor mil set to 

{-1.12013 - 0.261024 1, -0.492206 - 1.6121 1, 1.84975 + 1.77393 1, 2.39887 - 1.57675 1}. 
Momentum for spinor pl2 set to 

{-0.28478 + 0.261024 1, -0.713906 + 1.6121 1, -2.16198 - 1.77393 1, -1.74957 + 1.576 75 1). 

sol2hl = MHV[nai, pl4, 4, 5] MHV[1, ml2, pll] MHVbar[na3, pl2, 2] MHVbar[3, mli, pl3] // N // Chop 



sol2h2 = MHV[nill, pl4, 4, 5] MHVbar [nil2, pll, 1] MHV[pl2, 2, ml 3] MHVbar [3, nil4, pl3] // N // Chop 





We can verify that the result is as expected [87] 



is 1, 2 s 2, 3 

MHV[1, 2, 3, 4, 5] //N 

2 

-20.5335 + 31.3785 1 
sollhl + sollh2 + sol2hl + sol2h2 
2 

-20.5335 + 31.3785 1 



4-3 Spinor Integration 



In the following example we discuss how SOM can be used to perform the eval- 
uation of a double-cut by means of the spinor-integration method described 
in [42,48]. The Mathematica notebook file S@M_Spinor Integration. nb con- 
taining this example is distributed with the package. In particular, we consider 
the cut in the S34-channel of the one-loop four-gluon amplitude ^"^'^^^""^(l", 2+, 3~, 4+) 
with a gluon running around the loop. The S34-cut has two contributions, ac- 



40 



3- 

4+ 



Fig. 3. Double-cut for the determination of the coefficient of the bubble integral 
-^2('S34) in the one-loop 4-point gluon amplitude, with a gluon propagating in the 
loop. 

cording to the choice of the internal helicity. We will only take one of the two 
possibilities into account, namely 



where d$ is the standard Lorentz invariant two-body phase-space, see Figure 
3 We begin with the definition of the tree-level amplitude, and of the cut- 
integral 



GenMoinenta[{l, 2, 3, 4, 5, 6}] 




Momenta for the spinors 1, 2, 3, 4, 5, 6 generated. 




DeclareSpinor[ll, 12, A] 




{11, 12, A} added to the list of spinors 




DeclareSpinorMomentum[A, {Sqrt[3], 1, 1, 1}] 




Momentum for spinor A set to jVs", 1, 1, l|. 




DeclareLVectorMomentvrai[P34, Num4V[3] +Nvim4V[4]] 




{P34} added to the list of Lorentz vectors 




Four Momentum P34 set to 

{3.81711603574039234, 1.27297206351920934, -2.40158026409676221, 
mpmp[xl_, x2_, x3_, x4_] : = 

Spaa[xl, x3]''4/ {Spaa[xl, x2] *Spaa[x2, x3] *Spaa[x3, x4] *Spaa[x4, 


.39476806691977571} . 
xl] ) ; 



By using the momentum conservation £i — £2 = ^3 + ^4 = -P34, one can elimi- 
nate the dependence on £2, and write the integrand just in terms of £1 



Cut[l] = d$ * n53n53[l, 2, 11, 12] *inptnp[12, 11, 3, 4] 

dS <11 I 1)" (12 I 3)* 

<11 I 12)^ <11 I 2> <11 I 3> <12 I 1) <12 I 4) <1 I 2) <3 | 4> 
Cut[l] = Cut[l] /. { 

Spaa[x_, 12] -» Spab[x, P34, 11] /Spbb[12, 11], 
Spaa [12, x_] -» Spba[ll, P34, x] /Spbb[ll, 12] 

} 

- (dS <11 I 1>" <3 I P34 I 11]*) / 

(<11 I 2> <11 I 3) <1 I 2> <3 I 4) <11 I P34 | 11]^ <1 | P34 | 11] <4 | P34 | 11] ) 

Then we redefine the loop spinor variables according to the CSW [81] prescrip- 
tion: 




l^i) 



(38) 
(39) 
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J J J (AIP34IA] V (^1^^341 A] y ^ ' 

JdX^ J (AdA)[AdA], (41) 



|A]=|A>* 



Cut[l] = ASpinorReplace[Cut[l] , 11, Sqrt[t] *A] 

- (d56 <A I 1>^ <3 I P34 I 11]^) / (<A I 2> <A I 3> <1 I 2> <3 | 4) <A | P34 | 11]^ <1 | P34 | 11] <4 | P34 | 11] ) 
Cut[l] = BSpinorRfiplace[Cut[l] , 11, Sqrt[t]*A] 

- fd'!' ?. I 1^^'" '^3 I ?34 I ?.\^] ^ ' ?. I 2^' ■''A I 3'- ■''1 | 2^^ ^^3 | ' } \ ?34 | } ~ '1 \ ?34 | }.\ ''4 | ?3- | 3] 

We can perform the ^-integration by substituting the value of t as imposed by 

the cy-function . 

Cut[l] =Cut[l] //. dS-»dA*t/Spab[A, P34, A] 

- (dAt <A I 1>^ <3 I P34 I A]^) / (<A I 2> <A I 3> <1 I 2> <3 I 4> <A I P34 I A]^ <1 I P34 I A] <4 I P34 I A] ) 
Cut[l] = Cut[l] //. {t -» s[3, 4] /Spab[A, P34, A] } 

- (dAs3x4 <A I 1>^ <3 I P34 I A]*) / ((A I 2> <A I 3> <1 | 2> <3 | 4) <A | P34 | A]\l | P34 | A] <4 | P34 | A] ) 

Now we choose to integrate by parts in |A] and to extract residue in |A). To 
that aim we have to cast the integrand in a suitable form. We perform the 
reduction on the |A]-variable by means of the Schouten identities, 

Cut[l] = BSchouten[Cut[l] , A] //Expand; 

Cut[l] = Cut[l] //. (sp : (Spaa | Spab | Spbb) ) [x , P34, P34, y ] -»s[3, 4] * sp[x, y] 

dAs3x4 <A I if <A I 3> (1 I 3> dAs3x4 <A I 1> <A | 3> (1 | 3>^ 



<A I 2> <A I 4) <1 I 2) <A I P34 | A] <A | 2> <A | 4> <1 | 2> <3 | 4) <A | P34 | A]' 
dAs3x4 <A I 1>^ <A I 3> <3 I 4> dAs3x4 <A I 1>^ <A | 3> <3 | P34 | A] 



<A I 2> (A I 4) <1 I 2> <A I P34 | A] <A | 2> <A | 4> <1 | 2> <A | P34 | A] 

dAs3-4 <A I if (X I 3> <1 I 3> (3 | P34 | A] dAs3x4 <A | 1>^ <A | 3> <3 | P34 | A]^ 



<A I 2) <A I 4> <1 I 2> <3 I 4> <A | P34 | A]^ <A | 2> <A | 4> <1 | 2) <3 | 4) <A | P34 | A] ' 
dAs3x4 <A I 1> <1 I 3f dAs3x4 <1 I 3)^ 



(A I 2> <A I 4) <1 I 2> <A I P34 | A] <4 | P34 | A] <A | 2> <1 | 2> <3 | 4> <A | P34 | A] <4 | P34 | A] 
dAs3x4 <A I if <1 I 3> <3 I 4> dAs3x4 <A I 1>' <3 | 4>' 



<A I 2> <A I if (1 I 2> <A I P34 | A] <4 | P34 | A] <A | 2> <A | if (1 | 2> <A | P34 | A] <4 | P34 | A] 

dA S3 ^ 4 < 1 I 3 > * 



<A I 2) <A I 3> <1 I 2> <3 I 4> <1 | P34 | A] <4 | P34 | A] 
Cut[l, 4] = Coefficient [cut [1] , l/spab[A, P34, A] * l/ Spab [A, P34, A]* //Expand 

The integrand is thus expressed as sum of four terms according to their |A] 
dependence. 

[•A]" 1 

^ ^ n = 0, 1,2 or 



(A|P34|A]"+2 ' ' (A|P34|A](.|P34|A] 



42 



dAs3x4 <A I 1)^ <A I 3> <3 I P34 | A] ^ 

<A I 2> <A I 4> <1 I 2> <3 I 4> <A | P34 | A] ^ 

Cut[l, 3] = Coefficient [Cut [1] , l/Spab[A, P34, A] ^] * l/ Spab [A, P34, A]^ //Expand 

dAs3x4 <A I 1)^ <A I 3> <3 I P34 | X] dAsav, <A | (A | 3) <1 | 3> <3 | P34 | A] 
+ 

<A I 2> <A I 4)^ <1 I 2> <A I P34 | A] ^ <A | 2> <A | 4) <1 | 2> <3 | 4> (A | P34 | A]^ 
Cut[l, 2] = Coefficient [cut [1] , l/Spab[A, P34, A] ^] * l/ Spab[A, P34, A]^ //Expand 

dA S3x4 <A I if <A I 3> <1 I 3> 

<A I 2> <A I 4>^ <1 I 2> <A I P34 | A]^ 

dAs3v4 <A I 1> <A I 3) <1 I 3)^ dAs3 4 <A | 1>^ <A | 3) (3 | 4> 

<A I 2) <A I 4) <1 I 2) <3 I 4) (A | P34 | A]^ <A | 2) <A | 4)^ <1 | 2) <A | P34 | A] ^ 
Cut[l, 1] = Coefficient [Cut [1] , l/Spab[A, P34, A]] *1/ Spab[A, P34, A] //Expand 

dAs3x4 <A I 1> <1 I 3>^ dAs3v4 <1 I 3)^ 

h 1- 

<A I 2> <A I 4> <1 I 2> <A I P34 | A] <4 | P34 | A] (A | 2) <1 | 2) <3 | 4) <A | P34 | A] <4 | P34 | A] 

dAs3x4 <A I 1)^ <1 I 3) <3 I 4) dAs3x4 <A | 1)-^ <3 | 4)^ 

<A I 2> <A I 4>^ <1 I 2> <A I P34 I A] <4 | P34 | A] <A | 2> <A | 4>^ <1 | 2> <A | P34 | A] <4 | P34 | A] 
Trin4[l] = Expand[Cut [1, 4] //. dA->dea*deb] 

For brevity, we discuss only the integration of the term defined as Cut [1,4] . 

To keep track of the integration we use the exphcit definition of dX given 
above, dX = (A dX) [X dX] = dea * deb . After simphfying the integrand with 
trivial spinor identities, we integrate by parts in |A], using the identity 



^ ^ {X\Ps4\X]-+' (n+1) {X\Pu\X]-+HX\PM ^ ^ 

for n = 2 and |?7] = |4] , 

dea deb 33x4 <A I 1>^ <A | 3> <3 | P34 | A]^ 



<A I 2> <A I 4> <1 I 2> <3 I 4> <A | P34 | A] ^ 
IdySet = { 

Spab[4, P34, x_] ->Spaa[4, 3]*Spbb[3, x] , 

Spab[3, P34, x_] ->Spaa[3, 4]*Spbb[4, x] , 

Spab[x_, P34, 4] ->Spaa[x, 3] *Spbb[3, 4], 

Spab[x_, P34, 3] -»Spaa[x, 4] *Spbb[4, 3] 

}; 

Trm4[l] =Trm4[l] //. IdySet 

deadebS3x4 <A | 1>^ <A | 3> <3 | 4> [4 | A]^ 



<A I 2) <A I 4> <1 I 2> <A I P34 | A] 
IBPidy = 



debSpbb[4, A]^ Spbb[4, A]^ 



Spab[A, P34, A]^ 3Spab[A, P34, 4] Spab[A, P34, A]^ 

deb [4 I A] ^ [4 | A]^ 



<A I P34 I A]' 3 <A I P34 | A] ^ <A | P34 | 4] 
Tnn4[l] =Trm4[l] //. IBPidy 

dea 33x4 <A I l)-^ <A | 3> <3 | 4) [4 | A]^ 
3 <A I 2> <A I 4> <1 I 2) <A I P34 | A] ^ <A | P34 | 4] 



The final integration over |A) can be performed by summing over the residues 
at the simple poles in |A). The expression of Trm4[l] has two simple poles 
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at |A) = |2), |4). We notice that the residue at |A) = |4) vanishes, due to the 
term [4A]^ in the numerator. Therefore the result is given just by the residue 
at |A) = |2), 

Trin4[l] =Trm4[l] //. IdySet 

dea S3 4 <A I <3 | 4) [4 | A]^ 

3 <A I 2> <A I 4> <1 I 2> <A I P34 | X]-' [4 | 3] 

Two comments are in order. In this simple case, after the | A] -integration, the 
integrand contained only simple poles in |A). In general it may well happen 
that higher poles are present. Should this be the case, one can apply the func- 
tion ASchouten in order to single out the simple poles beneath the higher 

ones, and then take the residue in |A). 

The integration of Cut [1,3] and Cut [1,2] proceeds along the lines just out- 
lined. The result, as for Cut [1,4], will be a combination of rational functions 
of spinor products. Therefore they all contribute to the coefficient of the ln(s34) 
of the whole amplitude. 

The integration of Cut [1,1] is a bit different. In order to apply the integration 
by-parts in |A], one has to introduce a Feynman parameters. Then, the inte- 
gration over the spinor variables can be carried on as outlined above. Finally, 
one performs the integration over the Feynman parameter. The last paramet- 
ric integration is responsible for the rising of the logarithmic terms of the 
double-cut, whose coefficient can be directly assigned to the (poly-) logarithms 
of the whole amplitude. 



5 Conclusion and Outlook 

We have presented the Mathematica package SOM whose aim is to provide its 
user with a tool for performing the basic spinor algebra, plus the spinor-shifts 
needed for a very efficient analytic evaluation of scattering amplitudes at tree- 
and loop-level, accompanied by the support of the numerical evaluation at 
every computational stage. 

The basic properties of the functions introduced within SOM render it a very 
flexible program which could be further enriched with additional routines de- 
signed for more specific tasks. 
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A Functions Index 



The following table lists the functions provided by the package and the page 
where they are described. 



Command 


page 


Command 


page 


Command 


pa; 


$SpinorFunctions 


31 


Gaiifflia2 


31 


SmBA4 


20 


ACompactif y 


26 


GaimnaS 


31 


SmBA 


18 


ASchouten 


29 


Gamma5 


31 


Sm 


18 


ASpinorReplace 


29 


GenM omenta 


32 


SpOpen, SpClose 


23 


ASpinorShift 


30 


LVectorQ 


15 


Spaa 


20 


BCompactif y 


26 


Lat 


13 


Spab 


20 


BSchouten 


29 


La 


13 


Spba 


20 


BSpinorReplace 


29 


MP2 


15 


Spbb 


20 


BSpinorShift 


30 


MP 


15 


SpinorQ 


12 


CLat 


13 


Num4V 


33 


Sp 


12 


CLa 


13 


PfrinCSm2 


35 


To2DiinSpinor 


25 


CSm2 


19 


PfrmSin2 


35 


To4DiinSpinor 


25 


CSmBA2 


19 


PfrmSm4 


35 


ToSpinorLabel 


26 


Compactify 


26 


Proj Minus 


31 


USpa 


13 


Convert SpinorsToS 


23 


Pro j Plus 


31 


USpb 


13 


DeclareLVectorMomentum 


32 


s 


16 


UbarSpa 


13 


DeclareLVector 


14 


Schouten 


28 


UbarSpb 


13 


DeclareSpinorMomentum 


32 


ShiftBA 


30 


UnCompact 


27 


DeclareSpinor 


12 


SMatrixQ 


18 


UndeclareLVector 


14 


ExpandSToSpinors 


23 


Sm2 


19 


UndeclareSMatrix 


18 


GammaO 


31 


Sm4 


19 


UndeclareSpinor 


12 


Ganunal 


31 


SmBA2 


19 
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A list of all the functions is stored in the variable $SpinorsFunctions. 
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